Confinement Phenomenology in tiie Bethe-Salpeter Equation 
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We consider the solution of the Bethe-Salpeter equation in Euclidean metric for a qq vector meson 
in the circumstance where the dressed quark propagators have time-like complex conjugate mass 
poles. This approximates features encountered in recent QCD modeling via the Dyson-Schwinger 
equations; the absence of real mass poles simulates quark confinement. The analytic continuation in 
the total momentum necessary to reach the mass shell for a meson sufficiently heavier than 1 GeV 
leads to the quark poles being within the integration domain for two variables in the standard 
approach. Through Feynman integral techniques, we show how the analytic continuation can be 
implemented in a way suitable for a practical numerical solution. We show that the would-be qq 
width to the meson generated from one quark pole is exactly cancelled by the effect of the conjugate 
partner pole; the meson mass remains real and there is no spurious qq production threshold. The 
ladder kernel we employ is consistent with one-loop perturbative QCD and has a two-parameter 
infrared structure found to be successful in recent studies of the light SU(3) meson sector. 
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I. INTRODUCTION 

QCD models based on the Dyson-Schwinger equations 
[DSEs] provide an excellent tool in the study of non- 
perturbative aspects of hadron properties in QCD [Q. 
Such models can implement quark and gluon confinement 
[1,0,^, Q , generate dynamical chiral symmetry breaking 
[5,^, and maintain Poincare covariance. It is straight- 
forward to implement the correct one-loop renormaliza- 
tion group behavior of QCD 0, and obtain agreement 
with perturbation theory in the ultraviolet region. Pro- 
vided that the relevant Ward identities are preserved in 
the truncation of the DSEs, the corresponding currents 
are conserved. Axial current conservation induces the 
Goldstone nature of the pions and kaons § ; electromag- 
netic current conservation produces the correct hadronic 
charge without fine-tuning 

Previous work ^ has shown this model to provide 
an efficient description of the light-quark pseudoscalar 
and vector mesons via two infrared parameters within 
the rainbow truncation of the DSE for solution of the 
dressed quark propagators coupled to the ladder approx- 
imation for the Bethe-Salpeter equation [BSE]. Further- 
more, in impulse approximation, the elastic charge form 
factors of the pseudoscalars and the electroweak 
transition form factors of the pseudoscalars and vec- 
tors [|l2|, |l^ are in excellent agreement with data. As 
a coUorary, the strong decays of the vector mesons are 
well-described in impulse approximation without param- 
eter adjustment 0. 

Nonperturbative solutions of the DSEs are almost al- 
ways implemented in the Euclidean metric for technical 
reasons and this work is no exception. A resulting com- 
plication is that the solution of the BSE for meson bound 
states requires an analytic continuation in the meson mo- 
mentum to reach the on-mass-shell point = —AP. 
This causes the quark in the BSE to vary throughout 



a domain of the complex plane bounded by a parabola 
which widens in proportion to the mass M of the meson. 
For low mass mesons, such as the tt and K mesons, this 
excursion into the complex plane is simple to handle. Al- 
though the difhculties encountered in studies of ground 
state vector mesons are greater, they too may be over- 
come directly . For mesons with masses greater than 
about 1 GeV, these excursions into the complex plane 
are deep enough to encounter singularities in the rainbow 
DSE solutions for the quark propagators. As a result, lit- 
tle is known about the implementation of the BSE in this 
situation or even whether a solution is well-defined. For 
this reason, DSE-BSE model studies of mesons have for 
the most part been restricted to light mesons. An excep- 
tion is the study of Ref. |l5j in which the Ansatz adopted 
for the quark propagators in the non-Euclidean domain 
precluded singular behavior. 

It has long been known that the rainbow approxima- 
tion to the fermion DSE produces complex conjugate sin- 
gularities on the timelike half of the complex-momentum 
plane (3fJ(p^) < 0) in QED Thcphenomenon 
has been studied in detail within QED |l^ and in 
an Abelian-like model of QCD j20|. In these early stud- 
ies the singularities were attributed to artifacts of the 
rainbow approximation. A connection between the gen- 
eration of complex conjugate singularities and the pos- 
sible absence of real mass poles was noticed in a QED 
study Later studies within QCD |l] raised sug- 
gestions that complex singularities may not be artifacts 
of the approximation at all, but rather may be a prop- 
erty of the full theory, possibly related to confinement 
and the absence of real quark mass poles. Beyond the 
rainbow approximation, very little is known about DSE 
solutions and their singularity structure. These questions 
and the influence of approximations have yet to be com- 
pletely resolved. However they are being addressed; for 
example, in one simple model calculation, the inclusion 
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of the full dressing of the quark-gluon vertex in the DSE 
leads to the disappearence of all quark singularities in the 
complex plane save essential singularities at infinity ||]. 

It has been well-established by Dyson-Schwinger stud- 
ies that in order to obtain an acceptable amount of dy- 
namical chiral symmetry breaking in agreement with em- 
pirical observations, the effective rainbow DSE kernel 
must have a significant infrared strength. This in turn 
tends to move mass poles away from the real timelike 
axis, at least for a domain characterized by several 
GeV. This suggests that quark propagators may be suc- 
cessfully modeled in terms of parametrizations based on 
entire functions or at least based on functions with 
singularities away from the real axis such as complex- 
conjugate poles. 

If such an approach to modeling QCD is to be extended 
to include qq and qQ meson spectra beyond 1 or 2 GeV, 
then it is necessary to determine how solutions of the BSE 
are defined when complex singularities in the quark prop- 
agators may be encountered. This issue is addressed in 
the following using a simple parmetrization of the quark 
propagator. Although the method developed herein is 
first applied to the BSE, it has applications beyond me- 
son spectroscopy. For example, in production processes 
at moderate momentum transfers, the use of complex 
mass pole representations for confined quarks has been 
suggested as a means to handle the sizable timelike mo- 
mentum carried by the quarks in these processes [ p3| . 
One deals with these issues also in applications of the 
DSE-BSE approach to electromagnetic elastic and 
transition form factors |l^ of light mesons at high- 
momentum transfer. 

To simplify the analysis, we adopt a parameteriza- 
tion of the dressed quark propagator in terms of com- 
plex conjugate mass poles. In the standard approach 
to the meson BSE, two of the four integration variables 
are integrated numerically and it is during the integra- 
tion of these variables that poles may be encountered 
as a result of analytic continuation. In this article, a 
clear implementation, based on Feynman integral tech- 
niques, is introduced to map the four-dimensional inte- 
gration in the BSE onto a one-dimensional domain where 
numerical methods of evaluating Cauchy integrals are 
well-established. We find that when a complex conju- 
gate mass pole parametrization of the quark propagators 
is employed, the would-be qq decay width of the meson 
generated from a particular mass pole is exactly cancelled 
by the contribution of the conjugate partner pole. As 
a result, the meson mass M remains real valued, even 
above the pseudo-threshold > 3?(4m^) where mi is 
the lowest complex mass in the propagator representa- 
tion. 

The primary limitation of the method developed here 
is that it relies on an explicit parameterization of dressed 
quark propagators in terms of pairs of complex conju- 
gate mass poles. However the limitation is offset by 
the finding that such a representation can provide ex- 
cellent fits to both DSE solutions (as is shown in Sec. ^) 



and also to recent lattice-QCD data |24 . Of course, use 
of such a phenomenological representation will not ex- 
actly preserve the dynamical relation between vertices 
and propagators necessary to satisfy the Ward Taka- 
hashi identities of QCD. In particular, without the ax- 
ial vector Ward Takahashi identity, the chiral-limit pseu- 
doscalar meson solutions will not obey Goldstone's theo- 
rem [||. As a result, the toy model introduced in Sec. ^ 
would not be recommended for studies of light pseu- 
doscalar mesons; chiral-symmetry-preserving DSE ap- 
proaches have already been well-established |Q, |l^ for 
that sector. Rather, the present work explores techniques 
for extending the present DSE-based approach to mesons 
with masses above 1 GeV where one expects to encounter 
complex conjugate singularities. 

The article is organized as follows. In Sec. II, our 
notation for the ladder BSE for meson states is intro- 
duced, as is the complex mass pole representation of the 
quark propagators which will be employed. In Sec. Ill, 
a simplified BSE integral equation is used to illustrate 
the technique for analytic continuation of the BSE. The 
fact that the resulting bound state mass will be real is 
demonstrated explicitly. In Sec. IV, the method is ex- 
tended to the general BSE for a vector meson. A simple 
model for the ladder BSE kernel is introduced in Sec. V 
and numerical results are presented. A discussion of this 
approach and outlook is presented in Sec. VI. 



II. THE BETHE-SALPETER EQUATION 

In the following the Euclidean metric is employed. 
The scalar product of two four vectors is then 
a ■ b = and > for a spacelike 4- vector a^. 

The Dirac 7-matrices are Hermitian 7^ = 7^ and obey 
the anti-commutation relation {7^, 7^} = 25^jy. The 
dressed quark propagator S{p) and meson Bethe-Salpeter 
(BS) amplitude T{p; P) are solutions of the renormalized 
DSE, 



+Z, J g'D^,{p-q)-j^Siq)Tl{q,p),{l) 



and the BSE, 



.A 

r{p;P)^ / Kip,q;P)Q{Siq+)riq;P)Siq-)), 
J a 



(2) 



Here D^y{k) is the renormalized dressed-gluon propaga- 
tor, T\,{q,p) is the renormalized dressed quark-gluon ver- 
tex, r(p; P) is the BS amplitude for a quark- antiquark 
bound-state meson, and K{q,p;P) is the renormalized 
two-particle irreducible qq scattering kernel. We con- 
sider equal mass quarks with momenta q± — q ± P/2 
where q is the quark- antiquark relative momentum, and 
P = q^ — q^ is the meson momentum which satisfies 
p2 = ^Yieie M is the meson mass. 
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In Eq. (y|) the double contraction of Dirac indices has 
been denoted by 0. If the Dirac indices of the elements 
were reinstated in the BSE and written as parenthetic 
Roman letters (a, b, . . .), then Eq. ^ would appear as 

^iab){p;P) = / K(^ab;cd)iP,q;P)S(de){q+) 
Jq 

xr(,/)(q;P)%,)(g_) . (3) 

Use of the double contraction operator allows us to 
suppress Dirac indices, and simplify the forms of the 
equations appearing herein. 

In Eqs. (|^) and (^) = d'^k/ denotes a trans- 
lationally invariant ultraviolet regularization of the mo- 
mentum space integral with mass-scale A. Lorentz co- 
variance entails that the solution of Eq. (|^) is of the form 
= i^A{p^)+ Biyp^). The Lorenz scalars A{p^) and 
B{p^) are renormalized at spacelike — such that 
A{fi^) = 1 and B{^'^) = mq{^), where mq{fi) is the renor- 
malized current quark mass. After renormalization, one 
is free to remove the regularization by taking the limit 
A ^ oo. 

Recently, a successful model has been developed that 
provides an excellent description of the masses, decays 
and other properties of the light pseudoscalar and vector 
mesons [[zl, |ic|. The model consists of the rainbow-ladder 
truncations for the quark DSE (|^) and meson BSE (|^) 
and the use of an effective qq interaction constrained to 
coincide with perturbative QCD in the ultraviolet do- 
main and containing a phenomenological infrared behav- 
ior. In this particular model, the DSE kernel in rainbow 
truncation is 

Z,g^D^,ik)ri{q,p) ^ 4™eff(fc')i^;:r(fc)7.y , (4) 

where D^^f,°[k = p — q) is the free gluon propagator in 
Landau gauge, and acff(fc^) is the effective running cou- 
pling. The ladder truncation of the BSE kernel is 

K{p, q- P) -47raeff(fc2) Df^^k)^^^ , (5) 

where k — p — q and denotes the direct (or Carte- 
sian) product of matrices. This truncation scheme is 
self-consistent in that it ensures that the dressing of the 
quark-antiquark vector and axial-vector vertices gener- 
ated by Eqs. (^) and (||) satisfy their respective Ward- 
Takahashi identities. This feature is important since the 
axial- vector Ward-Takahashi identity guarantees that, in 
the chiral limit of zero current quark mass, the ground 
state pseudoscalar meson bound states are the massless 
Goldstone bosons arising from dynamical chiral symme- 
try breaking in accordance with the Goldstone's theo- 
rem p. The vector Ward-Takahashi identity guaran- 
tees electromagnetic current conservation for mesons and 
nucleons if the impulse approximation is used to describe 
the current in terms of quarks |ll[ . The rainbow-ladder 
truncation is particularly suitable for the flavor octets of 



light pseudoscalar and vector mesons where higher-order 
contributions to the quark-gluon skeleton-graph expan- 
sion have significant cancellations ||2^, Q and so may be 
safely neglected. 

With the ladder truncation from Eq. (^), the BSE for 
a light-quark vector meson is 

A(p2)r,(p;F) = 
.A 

/ X^,(fc)7^^(g+)rp(g;P)5(g_)7. , 

Jq 

(6) 

with k = p — q and 

Here T p is a transverse four vector. We will refer to the 
vector meson under study here as the p meson although, 
due to a number of simplifications made for illustration 
purposes, the particular model solutions considered here 
should not be taken as physical representations of the p 
meson. 

A linear eigenvalue A(P^) has been introduced into 
Eq. so that it will yield solutions over a continuous 
range of P^. A physical bound state of mass Afp corre- 
sponds to an eigenvalue of A(P^ — —^Q) — 1- To find 
such a solution, it is clear that Eq. ra) must be ana- 
lytically continued from Euclidean space where P^ > 
to Minkowski space where P^ < and the on-mass- 
shell condition P^ = — may be realized. In terms 
of the Euclidean four momentum P^, such an analytic 

continuation would correspond to P4 = i^P'^ + , 
from which it follows that quark momenta = q^ 
—Mp/4:±i^/q^MpZ that enter the BSE are complex 
functions of the integration variables q^ > and the 
direction cosine —1 < z < +1. Therefore, in order to 
carryout such an analytic continuation of the BSE (|^), 
the dressed quark propagators S{qj.) must be known in 
the complex-momentum plane within a parabolic region 
having the negative real (timelike) point —M^/A as the 
apex and extending symmetrically about the real axis. 

Studies within the rainbow-ladder truncation find 
that the non-analytic points of u- and d-quark propaga- 
tors that occur nearest to the origin q^ — are located 
at g2 = -0.207 ± i0.331 GeV^ while those of the s-quark 
propagator are found at q^ = -0.376 ± i0.602 GeV^ @. 
This finding also agrees with qualitative observations 
from an earlier study p7| . From the positions of these 
non-analytic points one can define a critical mass Mc = 
V— P^ > above which these singular points enter the 
domain of integration for the BSE. For the model of 
Ref. these critical masses are Mc — 1.09 GeV for uu 
mesons and Mc — 1.47 GeV for ss mesons. Clearly, the 
non-analytic points lie outside the domains of integration 
for ladder BSE calculations of ground state pseudoscalar 
and vector mesons, but may lie within the integration do- 
main for heavier mesons. This observation provides one 
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of the main reasons that many BSE studies of mesons 
focus primarily on light meson states. 

The appearence of poles in the domain of integration 
of the BSE raises two questions: 1) does the analytic con- 
tinuation of the BSE to timelike momenta < yield 
physically meaningful quantities, and if so, 2) by what 
methods does one obtain the solutions of this singular 
integral equation? The goal of this article is to explore 
these questions. To facilitate an initial investigation, it 
is convenient to adopt a parametrization of the quark 
propagator S{p) as a sum of N pairs of complex conju- 
gate mass poles. 



N 



Sip) = E 



n=l 



if - 



ij) - 



(8) 



where to„ are complex-valued mass scales and z„ are 
complex coefficients. Comparison with the general form 
Siji)'^ = ij)A[]?)^ B(j?), allows identification of the 
amplitudes A(p^), B(jp'\ and the quark running-mass 
M ) = B (p2 ) I A (j? ) . Equivalently, comparison with 
the general form S{p) — ~i;^avip'^)+ crs{p^) allows iden- 
tification of the scalar and vector amplitudes of the 
propagator, usip^) and ovip^) respectively. These are 
needed for the BSE. For example, the quark propagator 
of Eq. (P) has a scalar amplitude of the form 



N 



(9) 



and a vector amplitude crv(p^) which is obtained by re- 
moving each of the m„ and m* from the numerators. 

It is shown in Sec. ^ that this form of quark propaga- 
tor with = 3 provides an excellent fit to DSE solutions 
that employ realistic effective qq interactions. Further- 
more, this form has been used recently to provide a para- 
metric fit of the chiral-quark propagator obtained 
from lattice-QCD simulations psj P9|. 



III. ANALYTIC CONTINUATION: A SIMPLE 
EXAMPLE 

A general study of the analytic continuation of the BSE 
is quite involved and complicated by the Dirac algebra 
necessary to describe quark-antiquark states. It is pru- 
dent to first proceed by simplifying several of these com- 
plications and focus on the issue of analytic continuation 
of the eigenvalue A(P^) and the associated integrations 
over the complex-conjugate poles in the quark propaga- 
These restrictions are lifted in Sec. IV where the 



tor 

approach introduced here is extended to the general BSE. 

In the following, assume that the quark propagator 
is well-represented by a single pair of complex conju- 
gate poles (that is, = 1 in Eq. (|8|)). Furthermore, 
assume that the vector-meson BS amplitude may be 
represented by the single Dirac amplitude Tp{p]P) = 



(Spa — PpPa/P'^)jaV{q] P), whcrc Pp is the four mo- 
mentum of the bound-state vector meson. After these 
simplifications, the BSE (|^) may be converted to a scalar 
eigenvalue equation for the amplitude V{q; P) by multi- 
plication from the left by 7p followed by a trace over the 
Dirac indices to obtain 



\{P^)V{p-P)^ 



E 

a,6=±l 



fab{p,q,P) Viq-P) 



(10) 



Observe that the kernel of this equation is bilinear in 
the propagator amplitudes as{q±) and (Tv{q±) but only 
the denominators of these functions are displayed explic- 
ity in Eq. ([l^). The remaining factors from the quark 
propagators and effective qq interaction have been col- 
lected into the scalar function fab whose detailed form is 
unimportant to the present discussion but may be derived 
following the discussion of Sec. IV. In addition, the de- 



nominators have been written in terms of ^2 = uij^ _ 77^2 
and A2 = 2miimj, where m = mn + irrij is the complex 
mass scale parameter in the quark propagator in Eq. (|^). 
This form allows a parametrization of the square of the 
quark mass and its complex conjugate m^* in terms 
of the values a = ±1. 

In principle, the g-dependcnce of the integrand in 
Eq. ( p^ is only known once the solution V{q; P) has 
been obtained. For the present, consider the artificial 
situation in which the solution of the BSE is known to 
be of the form V — g (t){q'^) U{zq) where Zq = q-P/ \/ q^P^ 
and <l){q^) and U{z) are known functions. In this case, 
the constant g may be eliminated from the BSE and fol- 
lowing an integration over d'^p one obtains. 



(11) 



a,b=±l 



and 



Xab{P^) = 



d^q 



Fab{q,P) 



(2^)4 {ql + ;i2 + iaA2) (^2 ^^2^ ,5^2) 



(12) 

where Fab{q, P) is a non-singular, Lorcntz invariant func- 
tion of q and P. The objective is then to evaluate the 
four integrals Xab{P^) and show how they may be ana- 
lytical continued to negative values of P^ . This is carried 
out explicitly in the remainder of this section. 

When p2 ig continued to negative values, the quark 
momenta in Eq. ( p^ ) satisfy — (g^)*. The two fac- 
tors in the denominator of Xa.^a(P^) (that is, a = — 6 in 
Eq. ( |l2|) ) will thus be complex conjugates. As a result, 
for P^sufficiently timelike, there will always be some 
value of the integration variable q^ for which both fac- 
tors in the denominator of Xa,-a{P^) vanish simultane- 
ously. The result is a non- analytic point in A-| (p2) 

and A \~{P^)- Such behavior would usually be indenti- 

fied with the decay threshold for the vector meson into 
an asymptotically-free quark and antiquark pair. How- 
ever, in the following, it is shown that when the in- 
tegrals in Eq. (O) are added together the imaginary 
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parts cancel exactly thereby eliminating the possibility 
of quark- antiquark production thresholds. In contrast 
to this, the two factors in the denominator of Xaa{P^) 
never vanish simultaneously and consequently there are 
no corresponding non- analytic points in either A++(P^) 
or A__(P2). 

To analyze the behavior in detail, it is advantageous 
to restructure the integral in Eq. ( [I4 ) so that the ana- 
lytic continuation from the Euclidean domain > to 
the timelike domain < is straightforward. For val- 
ues of P^ > 0, one may combine the two factors in the 
denominator using the Feynman parametrization. 



1 

AB 



and obtain 



, 2 



d^q Fab{q,P) 



(27r)4 D{q,P,y)^ 



(13) 



(14) 



y=-l 




FIG. 1: The trajectory of pole position q^_{y) given by 
Eq. is shown as the dashed curve. The solid curve is 
the contour required to implement analytic continuation. 



where y = a — P, and 



p2 



D{q, P, y) - y{q.p+z — A') + iq' + — +fi'+i 



(15) 

One then eliminates the term linear in q in the denomi- 
nator by a shift of the origin q^ ^ qu — yPul'^- The re- 
sulting denominator depends on P^, g^, and y but not 
on q ■ P. After this shift, the function Fab{q,P) in the 
numerator still depends q ■ P. All other variables save q^ 
and y can be integrated out yielding. 



Xab{P') 



dy / dq 



q^F,,{q\y-P^) 
[q^~ql,{P\yW 



(16) 



where 



ql,{P\y) = -^{l-y^)^^?-'- {y{a - b) + a + b) A^, 

(17) 

is independent of g^, and 



Fab{q^y;P^ 
Here Zq if 



1 



+1 

dZq 

-1 



l-z^gFab{q-lP,P) 



(18) 

the direction cosine defined by 
The singular behavior due to 
the propagators is now a single second-order pole in the 
variable q^. Its location 9^6 (^^'2/) is given in Eq. ( |l^ ) 
and is not in the path of integration for spacelike P^ > 0; 
the integral is regular in Euclidean space. 

The y integration is performed after the integration 
over q^. Hence, for a given value of y, the ^^-integration 
domain is the positive real axis and the integrand has a 
pole in the complex plane whose location g^^(P^,?/) 
moves with P^. The location has a negative real part 
for P^ > and remains negative as P^ passes through 
zero and becomes moderately timelike. Only when 



P^ < — will the real part of the pole location be- 
come positive and have a chance of impinging upon the 
domain of the g^-integration. For that to happen, the 
imaginary part of ^^^(y) would have to vanish. As an- 
ticipated above, Eq. ( p7|) shows that this never happens 

when a — b\ for A++(P^) and A (P^) the two factors 

in the denominator of Eq. (^) are never zero simultane- 
ously. It is therefore trivial to carryout these integrations. 
However, for a = —b, the imaginary part of q'^i,{y) can 
become zero when the Feynman variable y = and care 
must be taken to properly handle such singular integrals. 

As a specific example, the trajectory of the pole loca- 
tion g^_(y) is shown in Fig. |^. It begins at (— ^^,A^) 
when y — —1, moves along the dashed curve as y in- 
creases, crossing the real axis when y = 0, and termi- 
nates at (— /i^,— A^) when y — 1. For y < 0, the cal- 
culation of the integral is straightforward. However for 
y > 0, the pole in q- has crossed the real axis and the an- 
alytic continuation of the integral A^ (P^) is defined by 

continuous deformation of the contour in a manner that 
ensures the pole does not cross the contour of integration. 
This is illustrated in Fig. |l|, for the deformed contour cor- 
responding to y « -1-1/2. Deformation of a integration 
contour in this manner is equivalent to adding the pole 
residue to the result of integration along the real q^ axis. 
If the residue term were neglected, one would obtain a 
discontinuous result due to the crossing of a branch cut 
onto the wrong sheet. 

For the integral A ^{P"^), the pole trajectory traces 

out the same path as in Fig. |l| but in the opposite di- 
rection; that is, ?/ > corresponds to the upper half q^- 
plane and the pole residue will be of opposite sign from 
A^ (P^)- After analytic continuation to P^ < 0, the re- 
sults for the four integrals can be expressed as 



Aafc(P' 



+ 1 



dy 



ds 



sFabis,y;P^) 
[■s-qlbiy^P')? 
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X / dy—{sFab{s,y;P^)) 



, (19) 



where Q{x) is the Heaviside step function. 

In the timchkc region, inspection of Eq. ( |l9| ) shows 
that individual integrals Xab{P^) are complex valued and 
satisfy 



(20) 



The imaginary parts cancel in pairs thereby producing a 
real result for the eigenvalue A(P^) = J2a b=± ^o,b{P^)- 
Notice that this does not mean that one may neglect 
the addition of residue terms required by the analytic 
continuation of the integrals. Rather, when > the 
residue (second) term of Eq. (|l^) is not purely imagi- 
nary and the sum of the real parts of the residues are 
required for the final result for X{P^)- Of course, these 
residue terms appear only for more timelike than the 
pseudo-threshold point P^ — — that is, for bound 
state masses Mp > 2/i. 

Before extension of this treatment to the general prob- 
lem, consider the connection between the above discus- 
sion of analytic continuation and the usual case encoun- 
tered above physical two-body thresholds. In particular, 
consider the limiting case of infinitesimal A for the sin- 
gle integral A (P^)- This corresponds to the standard 

case of employing elementary propagators with the con- 
ventional boundary conditions of a negative infinitesimal 
—ie added to the real mass. Above the two-body thresh- 
old P^ < —4^^, there would be qq production described 
by the residue term in A__(P^). It follows that a bound 
state calculation, based on such an integral alone, would 
result in a complex mass describing an unstable meson 
with a finite qq decay width. However, in the present case 
the presence of conjugate quark poles ensures that the an- 
alytic continuation introduces both positive and negative 
would-be widths that cancel exactly. The resulting BSE 
eigenvalue A(P^) in Eq. (10) will be purely real and there- 
fore bound state masses defined by A(P^ = —M^) = 1 
will be real and will describe a meson with no qq decay 
width. 



9m PF3+ -f^e^^avPlaquPp Fi 

qfi P5 + (T^„q^ Fq + (Ti_,^Py Fr 

qf^craf3qaPi3 Fs). (21) 



In a more concise form this is 



rp(9;P) 



8 

Q = l 



A"(9,P)P„(g;P) 



(22) 



where A" are the Dirac covariants in Eq. (|2l|); for ex- 
ample, A;J=i = {6pf, ~ PpPp/P'^) -fp. When all indices of 
a product of two Dirac covariants are contracted in the 
manner corresponding to the usual trace operation, one 
obtains 



trD[A^(g,P) A^(g,P)] - A«(g, P) A^(g, P), 



tap{q,P), 



(23) 



where t/^aiq, P) is a Lorentz scalar function of g^, P^ and 
q ■ P, and we have used the double-contraction product 
notation introduced in Eq. (||). 

Of course, one is free to choose any other set of Dirac 
covariants A"(q, P) as long as they are complete. In par- 
ticular, one is free to choose a set of covariants for which 
tap Saf3 in which case the resulting set of coupled in- 
tegral equations may be of a simpler form than those 
obtained herein. However, in practice the behavior of 
such covariants can be poorly behaved near q — 0, and 
can lead to numerical difhculties |0. Therefore, the set 
of covariants given by Eq. ( ^ ) is preferable and the lack 
of orthogonality is easily accommodated by tajs^q, P). 

In the previous section, it was shown that when the 
q^ and q ■ P dependence of the BS amplitude r^(q;P) 
is known, the analytic continuation of the BSE may be 
carried out in a straight-forward manner. In practice, 
the functional form of the BS amplitude is known only 
after the solution to the BSE has been obtained. The ap- 
proach outlined in the previous section may be general- 
ized by expansion of the functional dependence of the BS 
amplitudes Fa{q] P) in terms of complete sets of known 
functions of q^ and Zq = q - P/ ^JepP^ . Consider the ex- 
pansion 



IV. GENERALIZATION TO THE 
BETHE-SALPETER EQUATION 

The technique described in the previous section is gen- 
eralized to allow solutions to the homogeneous BSE for 
an arbitrary quark-antiquark bound state. To be spe- 
cific, we consider an equal quark-antiquark vector meson 
bound state. One can write the most-general BS ampli- 
tude for a massive vector meson as a sum of eight Dirac 
covariants and Lorentz-scalar amplitudes Pc(g;P), 



rp(9;P) 



p2 



{lf.Fi+qp iF2 



Fa{q;P) 



(24) 



where Ui{zq) are Chebyshev polynomials of the second 
kind satisfying the ortho-normality condition 

2 



dz Vr^^C/j(z) Uj{z) ^ 5,j , (25) 



and (j)s{q^) are a set of ortho-normal functions, that are 
specified later, and that satisfy 



dq' Mil Mq") = Sst . 



(26) 
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The required number of basis terms will be determined 
by convergence of the BSE solution with respect to the 
meson observables of interest. For example, one partic- 
ular choice of basis {(ps{q^)} which proves to be efficient 
for observables dominated by infrared physics may be 
inefficient for observables sensitive to the ultraviolet. In 
the present work, we take convergence of the meson mass 
as the criterion. 

For brevity, we define the operator Osi{p) which 
projects out the expansion coefficients g%{P'^) from the 
BS amplitude in Eq. (||). Using Eqs. (|2|) and (|2|), one 
finds 

gf.iP') = a,(p)<-i,b,P)A^'0r,(p;P) 

xt;i,(p,P)trD [Af{p,P) Tp{p;P)] .(27) 

The set of coefficients g%{P^) will be refered to as BS 
amplitudes since they are equivalent to the BS amplitude 
Tp{p; P) in that they contain all the dynamical informa- 
tion regarding the solution of the BSE. The expansion 
coefficients gfi{P'^) are functions of P^ only. The oper- 
ator Osi (p) is an integral operator that is introduced for 
brevity and notational convenience. 

When both sides of the BSE (^) are projected accord- 
ing to Eq. (Urt) one obtains 



{2n) 



K{p) 



QK(p,q-P)Q{S{q+)Tp{q-P)S{q-)) 



ABfStj 



^piP) 



QKip, q- P) (A^(q+)A^(g; P)A^(<Z-)) 

xa^(4)aB(d)0t('z')C/,(^g)4' (28) 

where the sums are over /3 e {f,...,8}, and i,j G 
{0, f , . . .}, and A, P e {V, S}. The second equality fol- 
lows by expanding the BS amplitude T p{q] P) according 
to Eqs. (^ ) and In Eq. ( psj) another set of Dirac 

covariants, A'^(g) ~ {—i'f-q, 1} with A — V,S, have 
been introduced to allow the propagator representation 

At this point, it is helpful to specify a coordi- 
nate system. We take Pp. = \/P^(0, 0, 0, 1), Pp — 

\/p^(0, 0, ^1 — Zp, Zp), and the integration momentum 

is represented as 



cos ( 
sin < 



l-z2 



(29) 



The integrand in Eq. ( Pq ) has no dependence upon angle 
0; the only dependence on direction cosine Xq is through 
the p ■ q dependence of the BS kernel K{p,q;P). In 
Eq. (^ ) we may therefore substitute 



d^q- 



dq^ q^ 



dz„ 



1 



dx„ 



(30) 



The integration over Xq in Eq. (|2^) produces a function of 
p'^,q^,Zp, Zq and P^ which can be expanded in the basis 
functions. That is, we define a quantity V by 



fdXq t-l{p,P) 



A" {p)QK{p,q;P) 



0(A^(<Z+)A^(g;P)A^(<Z_)) 



^<P.ip-'nizp) V:4f^(P^) Mq')UA^,) -(31) 

stij 



Upon substituting this back into Eq. 
tains a discretized form of the BSE, 



8h above, one ob- 



HP') g^^iP') = vr,f^r(^') 5:^r.AP') g^P'') , (32) 

where repeated indices are understood as being summed 
over. Here we have introduced the projection of the prod- 
uct of propagator amplitudes in the form 



dq^q^ I dz, 



.^\-zlcj,s{q^)U,{zq) 



^OA{ql)aB{ql) Uq^)U,{zq) . (33) 

Thus far, the BSE has been reduced to the discrete 
eigenvalue problem given in Eq. (^). Physical solutions 
are at A(P^ = ~Mp) = 1 from which one identifies the 
existence of a bound state vector meson of mass Mp. The 
only task remaining is the determination of the elements 
in Eq. (|3^ ) that make up the kernel of the eigenvalue 
problem. These are V and Q. The generalized interaction 
V is determined in terms of the kernel by the projection 
in Eq. 

Consider the explicit form of Q. Clearly there is a 
strong similarity between the form of Eq. ( |33| ) and the 
simple integral whose analytic continuation was consid- 
ered in Sec. [II. This is clarified by the expansion of the 



quark propagator amplitudes in terms of the mass pole 
terms from Eq. (|^) in the form 



CTAiq') 



N 

E E 

n— 1 a— ±1 



7Aa 



(34) 



where a=+i _ z„m„ for A ~ V or S, respec- 

tively, and Z^"^^^^ = z* or z^ml^, respectively. Then, 
Eq. (p3| ) can be written as 



5^AP') 



N 

E 



E 



A_ Mq^)U'd^,)Mq^)Ui{^,) 

(2^)4 (g2 + ^2 + ^aAl){ql + + ^feA^) 



(35) 
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These integrals in Eq. ( |3^ ) are of the type Aab(P^) al- 
ready explored in our discussion of analytic continuation 
in Sec. 



[II 



The only difference is that the number of 
integrals that have to be performed can be large since it 
depends on the number of basis terms for representation 
of the BS amplitudes via Eq. (^), and the number of 
conjugate pairs of mass poles for representation of the 
propagators. 



The same methods employed in Sec. Ill may be used 
here to calculate the integrals Gft^j{P)io^ all values 
of while accounting properly for the analytic con- 
tinuation. In particular, the Feynman integral method 
used here will lead to a denominator [q"^ — q^mabi^'^^ 2/)]^' 
where the pole location is 



Inmabiy) 



p2 



-ly- 



iM - bAl 



aAl +bAl 



(36) 



The integration over is carried out first, followed by 
the integration over the Feynman parameter y, and fi- 
nally the integration over Zg. Integrations over the two 
remaining angles in Eq. (pq) are trivial. In general, for 
a quark propagator parametrized in terms of N pairs of 
complex-conjugate poles, there are (2iV)^ pole trajecto- 
ries, half of which impinge on the integration domain 
and lead to residue additions. The imaginary parts of 



the integrals in Eq. (35) cancel in pairs as described in 
the previous section. 

In ladder approximation, the determination of V is 
straight-forward. There are never any poles encountered 
when is analytically continued to negative values. Of 
course, singularities may be encountered when contribu- 
tions beyond ladder are maintained in the BS kernel, and 
these must be treated in a way analogous to Q. This pos- 
sibility is not addressed in the present work. 



V. MODEL CALCULATIONS 

A. Ladder-Rainbow Kernel 

Here we use a particular model for the DSE-BSE ladde- 
rainbow kernel to numerically implement the preceeding 
developments. The model is specified by the "effective 
coupling" acff (fc^) and we employ the Ansatz [l^ 



47rQeff(fc^ 



47r2 D fc2 



iln 



1 



qcd; 



(37) 



with 7™ = 12/(33 - 2Nf) and Ti-s) = (1 - exp j^)/s. 

The ultraviolet behavior matches that of the 1-loop QCD 
running coupling Q!<,(fc^). The resulting solution of the 
ladder-rainbow DSE-BSE system of equations in the UV 



region generates the correct 1-loop perturbative QCD 
structure. The first term implements the strong infrared 
enhancement in the region < < 1 GeV^ phenomeno- 
logically required to produce a realistic value for 
the chiral condensate. With rrit = 0.5 GeV, t = e^ — 1, 
Nf = 4, Aqcd = 0.234 GeV, and a renormalization scale 
/i = 19 GeV, it has been found that uj — 0.4 GeV and 
D = 0.93 GeV^ give a good description of (qq), m^/x and 
fjr with physically acceptable current quark masses, 
[lO| . The propagator amplitudes from this rainbow DSE 
model have recently been shown to have the same 
qualitative features as lattice QCD simulations |^8[ pgf . 

The vector meson masses and electroweak decay con- 
stants produced by this model are in good agreement 
with experiments pO| . Without any readjustment of the 
parameters, this model agrees remarkably well with the 
most recent Jlab data [ plj for the pion charge form fac- 
tor F7r((5^). Also the kaon charge radii and electromag- 
netic form factors are well described |l^, ^ . The strong 
decays of the vector mesons into a pair of pseudoscalar 
mesons are also well-described within this model [|l^, ^ . 



B. Results 

The DSE solution for the u/d quark propagator from 
the model interaction given in Eq. (^7|) can be well fit 
with = 3 pairs of complex conjugate poles in the repre- 
sentation given in Eq. (||). The result displayed in Fig. ^ 
corresponds to the parameter set 



mi = 0.547-l-i0.303 GeV, 

zi ^ 0.200-1-^0.475, 

ma = -1.262 -hiO.570 GeV, 

Z2 = 0.142 + i0.045, 

mg = 1.560 + ^0.564 GeV, 

Z3 = 0.160 + i0.015 . 



(38) 



This parameterization was chosen so that Z(p'^) and 
M{p'^) were reproduced well at both = and — 
100 GeV^ while the main features of the momentum 
dependence are preserved. We imposed the constraint 
that M(p^) should approach its UV limit from above. 
The pole locations are p^ = —mf where mf = ^"^ + iAf. 
Thus we have m = 0.455 GeV, Ai = 0.576 GeV, = 
1.13 GeV, A2 = 1.12 GeV, and = 1-46 GeV, A3 = 
1.33 GeV. 

With use of this propagator representation in the BSE, 
a, u/d quark meson bound state with mass greater than 
the lowest pseudo-threshold 2/ii =0.910 GeV would be 
needed to test the method we have described. The 
model-exact p meson mass from the kernel Eq. (|37| ) is 
0.742 GeV poj and is therefore not suitable. Heavier 
u/d mesons have not been carefully studied within this 
model interaction due to the very singularity issue we are 
concerned with. In fact the evidence from exploratory 
studies is that, for example, m^^ in this ladder-rainbow 
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model is only about 200 MeV above nip and also does 
not provide a clear test. 
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FIG. 2: The quark propagator running mass M{(f) and 
wave-function renormalization function Z{q^). The DSE cal- 
culation (circles) is compared to the fit in terms of three 
pairs of complex conjugate mass poles given in Eq. (p^. The 
toy model used to test the Bethe-Salpeter method uses two 
pairs of poles given in Eq. (Eol). 



To test the present method, we artificially modify the 
pole parameterization of the quark propagator so that the 
lowest pseudo-threshold moves below the vector meson 
mass, and the latter moves to 0.770 GeV. We achieve this 
with the following parameterization in terms of = 2 
pairs of conjugate mass poles 

mi = -0.400 -f- i0.200 GeV, 

zi = 0.15, 

m2 = 0.550 + ^0.350 GeV, 

Z2 = 0.35 + i0.37 . (39) 

The resulting propagator amplitudes are also shown 
in Fig. Ij. This parameterization corresponds to /ii = 
0.346 GcV, Ai = 0.4 GeV, and ^2 = 0.424 GeV, A2 = 
0.621 GeV. With this propagator, the lowest pseudo- 
threshold in the BSE is at 2^i =0.692 GeV. The sub- 
sequent BSE solution is not to be taken as a physical 
represenation of the p meson; its purpose is to provide 
a test calculation where the propagator singularities are 
within the domain of integation. 

The basis functions used to expand the dependence 
of the BS amplitudes are 



A 



Ls[w) e 



(40) 



where w = q^/h?', the typical hadron scale is A = 
0.800 GeV, Ls{w) are Laguere polynomials, and Ng is 




P" (GeV ) 

FIG. 3: The eigenvalue A(P^) for the vector meson BSE. 
The mass-shell at A(— rrip) — 1, indicated by the vertical line, 
is above the lowest pseudothreshold set by the complex mass 
poles of the quark propagators. Proper analytic continuation 
of the integral produces the solid line; analytic continuation 
of just the integrand produces the dashed line. 



fixed by the normalization condition given in Eq. (|26|). To 
obtain better than 1% accuracy for the eigenvalue A(P^), 
we found that 7 terms in this basis were sufficient. The 
number of Chebyshev terms used for the angle basis was 
one. The first covariant from the set given in Eq. (|2l|) is 



known to be dominant for the vector meson solution |1C 
and this was the only one used here. 

The resulting BSE eigenvalue A(P^) is given in Fig. ^ 
as the solid curve. In the spacelike region {P^ > 0) 
the eigenvalue is a monotonically, decreasing function 
of P^. When P^ decreases below zero and becomes 
timelike, the eigenvalue A(P^) increases for a while then 
crosses unity, at P^ = — M^, signifiying a bound state 
at Mp = 0.770 GeV. In Fig. | this value of jg shown 
as a thin, vertical line. On the timelike side of this line, 
the eigenvalue continues to increase, then begins to de- 
crease and ultimately goes to zero for (P^ ^ 0). This 
is because the eigenvalue is a measure of the magnitude 
of the BSE kernel and the propagators therein falloff as 
the difference between the momentum arguments and the 
pole positions becomes larger than any mass scale in the 
interaction. 

The solid curve in Fig. ^ represents the proper analytic 
continuation of the BSE from P^ > to P^ < 0. The 
dashed line represents the result from a direct evaluation 
of the BSE integral without regard to the possibility of 
singularities and their movement as a function of P^ in 
relation to the integration domain. That is, only a tran- 
scription from P^ > to P'^ < has been made in the 
integrand of Eq. (|35|); the Feynman technique for com- 
bining the denominators and shifting the 4-momentum 
variable to complete the square has not been made. The 
integration variables are q^ and Zq = q ■ P. It is clear 
that for the spacelike domain (P^ > 0), and for the lim- 
ited timclike domain for which no poles move into the 
domain of integration (— 4/if < P^ < 0, i.e., below the 
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lowest pseudo-threshold), the two approaches give iden- 
tical results. The dashed line represents the approach to 
the BSE that has usually been implemented; it has been 
limited to states light enough to avoid the propagator 
singularities. 

With the parameters given above, the quark poles 
are first encountered within the BSE integration when 
the analytic continuation reaches = — that is, 
at p2 ^ -0.48 GeV^ Although the contibutions of the 
conjugate pair of poles to the imaginary part of the BS 
eigenvalue A(P^) cancel, their real contribution produces 
a discontinuity in the derivative A'(P^) at that point. 
Careful inspection of the solid curve in Fig. ^ reveals 
that this cusp occurs at the precise value where the full 
calculation and the nai've calculation diverge. 

VI. DISCUSSION 

We have provided a method to obtain meson bound 
state solutions from the Bethe-Salpeter equation in Eu- 
clidean metric when the dressed quark propagators have 
time-like complex conjugate mass poles within the in- 
tegration domain. This approximates features encoun- 
tered in recent QCD modeling via the Dyson-Schwinger 
equations; the absence of real mass poles simulates quark 
confinement. For this exploratory study we represent the 
quark propagators as a sum of complex-conjugate mass 
poles, and project the BSE on to complete basis sets 
to represent the momentum and angle dependence. We 
use Feynman integral techniques to combine the propa- 
gator denominators and map the integration to a one- 
dimensional domain. This allows a clear analysis of the 
analytic continuation in total meson momentum needed 
to reach the mass shell and which causes the singulari- 
ties to impinge upon the integration domain. The BSE 
linear eigenvalue remains real; in other words the eigen- 
mass remains real. The would-be decay width from one 
pole is exactly cancelled by the effect of the partner pole; 
the meson is stable against decay into a quark-antiquark 
pair. This describes the confinement of quark and an- 
tiquark within the meson bound state, even though the 



meson mass is "above threshold" . 

One of the limitations of the method presented here 
is that it relies upon the projection of the BSE on to a 
complete basis set of functions to represent the momen- 
tum dependence of the BS amplitudes. In general the 
behavior of the BS amplitudes is not known until after 
the BSE is solved. If the chosen basis set is an ineffi- 
cient representation, one would expect this to result in 
a lack of convergence. The (vector) meson mass that we 
have studied here is an integrated quantity dominated by 
infrared physics; it is not surprizing that convergence is 
easily achieved with the basis states of Eq. ( ^0| ) which 
fall off exponentially in the ultraviolet region. However, 
we know from studies such as Ref. |l^ that the leading 
ultraviolet power law behavior of the vector meson BS 
amplitudes is l/g^. Clearly the present basis would be 
inefficient for applications that depend strongly on the 
BS amplitude in this domain, such as the asymptotic 
behavior of the pion charge form factor For such 

studies, appropriate basis sets would have to be utilized. 

The principal reason that the BS amplitudes are ex- 
panded in terms of known basis functions is that the 
Feynman integral technique for handling the two propa- 
gator denominators requires a complex variable shift for 
the rest of the integrand. In principle, the analytic con- 
tinuation of the BS amplitude into the complex plane is 
determined by the dynamics of the BSE and can be deter- 
mined only through the solution. The implementation of 
the present method of solution requires a known analytic 
behavior for the BS amplitudes as expressed through the 
basis functions. One would expect an inadequacy of the 
basis in this respect to show up as poor convergence. 

Future work using the BSE solution method presented 
here will include extension of the DSE modeling approach 
to SU(3) flavor meson states above 1 GeV and meson 
form factors for > 3 GeV^ 
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